21:46 Liquid Crystals manuscript 



Liquid Crystals 

Vol. 00, No. 00, Month 200x, 1-16 
INVITED ARTICLE 

Shape- dynamic growth, structure, and elasticity of homogeneously 
oriented spherulites in an isotropic/smectic-A mesophase transition 

Nasser Mohieddin Abukhdeir* and Alejandro D. Rey 

Department of Chemical Engineering, McGill University, 3610 University St., Montreal, 
Quebec, Canada 



A Landau-de Gennes model that integrates the nematic quadrupolar tensor order parameter 
and complex smectic-A order parameters is used to simulate the two-dimensional growth of 
an initially homogeneous smectic-A spherulite in an isotropic matrix. These simulations are 
performed in the shape-dynamic (nano-scale) regime of growth under two material condi- 
tions: isotropic nematic elasticity and equal splay-bend nematic elasticity. A comparison of 
the growth kinetics, spherulite morphology, interfacial/bulk energy landscapes between both 
cases is made showing that equal nematic splay-bend elasticity is required to reproduce past 
experimental and theoretical observations. Additionally, a previously unknown undulation 
instability during spherulite growth is found which, in conjunction with preferred planar an- 
choring and defect shedding mechanisms at micron length scales, could explain the formation 
mechanism of focal conic curvature defects and ultimately smectic-A "batonnet" structures 
observed experimentally. 
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1. Introduction 



The study of liquid crystalline materials has had a rapid and profound technological 
impact on civilisation in the past century. In the overarching field of Soft Matter, 
the study of liquid crystals has contributed fundamental advances in technological 
areas such as display technology and high-performance materials. Additionally, 
interest is quickly emerging in an area in which the understanding of liquid crystals 
is being recognised as crucial: biological systems (lil), and Nature as a whole. Self- 
organisation is a basis on which all biological systems have developed and evolved 
and the presence of liquid crystal phases in the very building blocks of Nature, the 
cell's phospholipid bilayer, is self-evident (0, d). As our understanding of biological 
systems increases .more and more instances of liquid crystal phases are found Q): 



in muscle tissue (U), the development of bone fljj), and even in the beginnings of 
life itself (0). 

Over the course of almost four decades, the vast contributions of P.G. de Gennes 
to the area of liquid crystals have enabled the advancement of knowledge up to 
this point. One of these contributions was his theoretical work on the smectic-A 
mesophase. His discovery, simultaneously with McMillan, of the analogy between 
superconductors and smectic-A phase ordering has opened up an avenue for study 
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of this mesophase via modelling and simulation. The key role that this mesophase 
plays in biological systems alone shows the importance of this single contribution 
of de Gennes. In his Nobel lecture in 1991, he states that "smectics bring me 
naturally to another important feature of complex fluids - namely that, in our 
days, it is sometimes possible to create new forms of matter" which alludes to the 
types of applications of smectics that could be developed in the future and to their 
roles in biological systems. 

The most simple of the smectics is the smectic-A mesophase, which displays 
lamellar translational ordering, in addition to the orientational ordering of nemat- 
ics. Recently an increasing amount of interest in this mesophase, in particular of 
materials exhibiting a direct isotropic/smectic-A (disordered/ordered) transition, 
has resulted in many experimental and theoretical results. Nonetheless, the under- 
standing of this mesophase is in a nascent stage. Much of this is due to the time 
and length scales at which the structures and dynamics occur being on the nano- 
scale. These properties provide a obvious application of theoretical study through 
modelling and simulation in order to both enhance experimental research and make 
predictions independently. 

A fascinating range of liquid crystal growth morphologies have been the focus 
on much study The smectic-A mesophase, with its lamellar ordering on the 
molecular scale, exhibits growth, defect, and texture phenomena not seen in the 
nematic phases, including the lamellar-like cholesteric mesophase. Focusing on the 
growth phenomena solely, on transition from the isotropic/disordered phase, a va- 
riety of self-assembled smectic-A structures have been observed to form. These 
unique morphologies can be attributed to complex dynamics involving interfacial 
tension anisotropy, which results in preferred anchoring, and bulk texturing. 

The current study is a part of an overall effort to understand kinetics, dy- 
namics, and morphology of the direct isotropic/smectic-A liquid crystalline phase 
transition. The theoretical focal point of this research is a high-order Landau-de 
Gennes type phenomenological model of Mukherjee, Pleiner, and Brand U Q)- 
This high-order model incorporates much of the key physics involved in the direct 
isotropic/smectic-A transition which occur on multiple scales and involve multiple 
types of phase-ordering: orientational and translational. In this effort a compre- 
hensive approach has been developed for the determination of phenomenological 
parameters for the model and efficient phase diagram computation Following 
this, numerical simulation was used to study phase transition kinetics and defect 
dynamics (lid ), surface effects Aim , and most recently spherulite growth (1 M& ). In 
addition to these simulation studies, theoretical contributions have been made to 
the study of smectic-A filamentary growth and buckling (173 . Q). 

Past work on smectic-A spherulite growth has focused on initially radially ori- 
ented nuclei (I MS ). The current work studies an alternate initial configuration where 
the nucleus is homogeneously, or "ideally," oriented. This and past work 
focuses on rod-like low molecular mass liquid crystals and phenomenological param- 
eters are based, in part, on experimental data from 12CB (dodecyl-cyanobiphenyl) . 
This work neglects nucleation mechanisms, thermal fluctuations, heat of transition, 
impurities, and convective flow while taking into account energetically the inter- 
coupling between orientational/translational order and variation of smectic layer 
spacing. The objectives of this simulation study are: 

(1) Determine the dynamic growth morphology of an initially homogeneous 
smectic spherulite with no preferred anchoring at the isotropic/smectic-A in- 
terface (isotropic nematic elasticity). 

(2) Determine the dynamic growth morphology of an initially homogeneous 
smectic spherulite with preferred planar anchoring (as is observed experimen- 
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Figure 1. schematics of the (a) isotropic, (b) nematic, and (c) smectic-A phases. 



tally) at the isotropic/smectic-A interface (equal bend-splay nematic elastic- 

ity). 

(3) Compare these simulation predictions with simplified shape equation mod- 
els. 

The paper begins with brief background on liquid crystals and phase-ordering 
transitions. The model is then presented and details of the simulation method/ 
conditions detailed. A brief discussion of the connection between two- and three- 
dimensional simulation is then given using a less complex example of initially ho- 
mogeneous spherulite growth in the isotropic/nematic transition. Finally, results 
of the work are presented and conclusions made. 



2. Background 

2.1. The first- order isotropic/smectic-A transition 

This work focuses the study of rod-like thermotropic liquid crystals which exhibit 
a first-order isotropic/smectic-A mesophase transition. An unordered liquid, where 
there is neither orientational nor translational order (apart from an average inter- 
molecular separation distance) of the molecules, is referred to as isotropic. Liquid 
crystalline order involves partial orientational order (nematics) and, additionally, 
partial translational order (smectics and columnar mesophases). The simplest of 
the smectics is the smectic-A mesophase, which exhibits one-dimensional transla- 
tional order in the direction of the preferred molecular orientational axis, which 
can be thought of as layers of two-dimensional fluids stacked upon each other. 
Other more complex types of smectics exist, for example tilted smectic/smectic-C 
and hexatic smectic/smectic-B mesophases. In this context, the relative simplicity 
of the smectic-A mesophase makes it the ideal starting point and, subsequently, a 
template phase for the vast set of self-assembled lamellar systems (15). Schematic 
representations of these different types of ordering are shown in figure [TJ 

Due to the first-order nature of the isotropic/smectic-A transition, a coexistence 
temperature interval exists where both the isotropic and smectic-A phase are either 
stable or metastable. The phase diagram computation method for the model used 
in this work was developed previously (5) and the resulting phase diagram for the 
phenomenological system used in this work is presented in figure El 
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Figure 2. Computed phase diagram for 12CB (0) with nematic scalar order parameter (solid line), smectic 
scalar order parameter (dotted line), and the normalised wave vector (stippled line). The isotropic/smectic- 
A coexistence region is enclosed by the vertical dashed lines, indicating the temperature range over which 
the isotropic/smectic-A phases coexist in stable/meta-stable states ([3. \Jdti. 



3. Modelling and Simulation 

3.1. The Landau- de Gennes model 

Theoretical characterisation of mesophase order is accomplished using order pa- 
rameters that adequately capture the physics involved. These order parameters 
typically have an amplitude and phase associated with them. In order to charac- 
terise the partial orientational order of the nematic phase, a second order symmetric 
traceless tensor can be used (Q): 

Q = S (^nn - -I^j + ^P (mm - 11) (1) 

where n/m/1 are the eigenvectors of Q, which characterise the average molecular 
orientational axes, and S/P are scalars which characterise the extent to which the 
molecules conform to the average orientational axes (0, 0, !£)• Uniaxial order 



is characterised by S and n, which correspond to the maximum eigenvalue (and 
its corresponding eigenvector) of Q, S = |/i n . Biaxial order is characterised by 
P and m/Z, which correspond to the lesser eigenvalues and eigenvectors, P = 

-| (Mm M/)- 

The smectic-A mesophase has one-dimensional translational order in addition to 
the orientational order found in nematics. Characterising this mesophase can be 
accomplished through the use of primary (orientational) and secondary (transla- 
tional) order parameters together (ldi ). A complex order parameter can be used to 
characterise translational order Q): 

^ = ^ (2) 

where <j) is the phase, ij) is the scalar amplitude of the density modulation. The 
density wave vector, which describes the average orientation of the smectic-A den- 
sity modulation, is defined as a = V<j)/\V<j)\. The smectic scalar order parameter 
i\) characterises the magnitude of the density modulation, and is used in a dimen- 
sionless form in this work. In the smectic-A mesophase the preferred orientation of 
the wave vector is parallel to the average molecular orientational axis, n. 

A Landau-de Gennes type model for the first order isotropic/smectic-A phase 
transition is used that was initially presented by Mukherjee, Pleiner, and Brand 
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|3 ; B) and later extended by adding nematic elastic terms (~2(\^2^): 



f-h = \a{Q:Q)-\b(Q-Q):Q + \c(Q: Q) 2 + \a |*| 2 + |*| 4 

-\W 2 (Q ■ Q) ~ \eQ : (V*) (V**) 

+1j 1 (VQ;VQ) + ^h(V Q V Q) + : (VQ : VQ) 

+ ^i|V*| 2 + ^ 2 |V 2 *| 2 
2 4 1 



(3) 



a a (T - T NI ) 
a a (T - Tai) 

where / is the free energy density, fo is the free energy density of the isotropic phase, 
terms 1-5 are the bulk contributions to the free energy, terms 6-7 are couplings 
of nematic and smectic order; both the bulk order and coupling of the nematic 
director and smectic density- wave vector, respectively. Terms 8-10/11-12 are the 
nematic/smectic elastic contributions to the free energy. T is temperature, Tni/Tai 
are the hypothetical second order transition temperatures for isotropic/nematic 
and isotropic/smectic-A mesophase transitions (refer to (ldi ) for more detail), and 
the remaining constants are phenomenological parameters. 

The Landau-Ginzburg time-dependent formulation (0) is used to capture the 
dynamics of the phase transition. Due to the higher order derivative term in the 
free energy functional, a higher order functional derivative must be used. Addition- 
ally, in order to utilise standard numerical solution techniques, the complex order 
parameter equation ([2]) is separated into its real and imaginary contributions (0): 

* = A + Bi (4) 

The general form of the time-dependent formulation is as follows (0): 




where /i r //i s is the rotational/smectic viscosity, and V the volume. As previously 
mentioned, a higher order functional derivative must be used due to the second- 
derivative term in the free energy equation ([3j): 

SF _df d ( df \ d d ( df \ 

where 9 corresponds to the order parameter. 

Substituting equation (j4j), the free energy ([3j), and high order functional deriva- 
tive (j7|) into the time-dependent formulation ([5]) yields the closed set of model 
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equations: 



dQ 

dt 



¥ Q -b*(Q- Q) ST + c*(Q:Q)Q-8* (A 2 + B 2 ) Q - -e* ( VAVA + VBVB 



+V • (ZJVQ) 

[a* A + (3* (A 2 + B 2 )A- S* A (Q : Q)] + V • 6J VA - e*Q ■ VA - -b* 2 V (V 2 A) 



.dA 

dt ~~ 

= - [a*B + (3* (A 2 B 2 ) B - S*B (Q : Q)] + V 



b\VB - e*Q ■ VB - -b* 2 V (V 2 B) 



(8) 



where the asterisk denotes an nondimensionalized value, the superscript ST denotes 
the symmetric/traceless portion of a tensor, and /x* is the ratio of the smectic and 
rotational viscosities. The nondimensionalized model parameters are as follows: 
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(9) 



where / is the simulation-specific imposed length scale. 



3.2. Simulation method and conditions 

A square computational domain with an imposed length scale of / = 2.93 x 10 _1 /xm 
(approximately 75 smectic layers at 330K (25), see figure EJ) was used in two sepa- 
rate simulations with isotropic nematic elasticity (no preferred interfacial anchor- 
ing Zi > 0, h,h — 0) and equal bend/splay nematic elasticity (preferred planar 
interfacial anchoring Zi, Z2 > 0, Z3 = 0). Referring to figure El Neumann boundary 
conditions were used to simulate bulk conditions. The initial condition for both 
simulations was a smectic-A spherulite in an initially homogeneous layer configu- 
ration (see figure 03). The radius of the spherulite was initially set to ro = 4.0nm. 
The initial value used for S, -0, and the layer spacing correspond to the homoge- 
neous values at T = 330K, determined from the computed phase diagram (figure 
ED- The Heaviside step function was used to generate the initial spherulite. The 
constraint that the spherulite does not impinge on the domain boundaries was 
verified post-simulation. 

A commercial package, Comsol Multiphysics, was used to solve the time- 
dependent model Quadratic Lagrange basis functions were used for the Q- 
tensor variables and quartic Hermite basis functions used for the complex order 
parameter components. Standard numerical techniques were utilised to ensure con- 
vergence and stability of the solution. This platform does not support adaptive 
mesh refinement, thus a uniform mesh was used with a density of approximately 
14.8 nodes/nm 2 . Previous simulations using this model and numerical method 
have shown good agreement with both past experimental and theoretical find- 
ings (fjfil . [Til ) . Additionally, exhaustive past work using this numerical method and 
the Landau-de Gennes model for the first-order isotropic/nematic phase transition 



February 9, 2009 21:46 Liquid Crystals manuscript 



Figure 3. Surface plot of Re(^) over the computational domain with the initial spherulite condition 
indicated by stippled circular region (assumed homogeneously oriented smectic-A nucleus). White/black 
corresponds to the maximum/minimum of Re(ty) and the imposed length scale is / = 2.93 x 10 _1 /im. The 

12CB, are T NI = 322. 85K, T AI = 
a = 1.903 x 10 6 — i--, (3 = 

1 x K)- 12 pN, 6i 



material parameters and phenomenological coefficients, based upon 
330. 5K, a = 2 x 10 5 — b = 2.823 x 10 7 ^4-, c = 1.972 x 10' 



3.956 x 10 



8 J 



, s 



m 3 K ' 

9.792 x 10 6 ^4> 



1.938 x lO- 11 ^, h 



1 x 10~ 12 pN, 



b 2 



3.334 x lO -30 Jm, fi N 
25. 



,4 x 10 



-2 Nxs 



, and the ratio of the rotational and diffusional viscosities 



used was 



tt2ui42yi) has served to further validate this simulation approach. 



3.3. Multidimensional computation 

Two fundamental challenges to research utilising numerical simulation are compu- 
tational limitations and functionality of numerical routines. As a result of the first 
challenge, the vast majority of simulation studies in the field of liquid crystals et 
al have been limited to one- and two-dimensions. This is the case for the current 
study as well, but the symmetries and similarities between full three-dimensional 
simulation and two-dimensional simulation (of three dimensional phenomena) pro- 
vide a strong motivation for and utility of obtaining these two-dimensional so- 
lutions. A specific justification for this relationship can be found by comparing 
past simulation results of a growing nematic spherulite in two-dimensions 
\2S ) with results of the same system (model, parameters, and initial conditions) in 
three-dimensions. Figured! shows results of three-dimensional simulation of the the 
growth of an initially homogeneous nematic spherulite based on two-dimensional 
studies by Wincure and Rey (I^6l4^). Figures S^-c show three orthogonal views of 
an initially homogeneous nematic spherulite where the spherulite morphology in 
the y/z- plane (figure UK) shows qualitative similarities to past two-dimensional re- 
sults (fliL \2(h. Qualitatively, both the two-dimensional and three-dimensional initial 
conditions share symmetries in one orthogonal plane, thus the two-dimensional so- 
lution is qualitatively a subset of the three-dimensional result. Based on this result, 
the two-dimensional smectic-A spherulite computations presented later are also ex- 
pected to contain signification features found in three-dimensional structures. 



4. Results and discussion 

Two different general types of post-nucleation growth are observed in growth pro- 
cesses of liquid crystals. The shape-dynamic regime of growth from an initial nu- 
cleus, dominated by bulk elastic energy, involves shape and texturing dynamics 
as interfacial anchoring effects become important when the spherulite radius sur- 
passes the characteristic length of liquid crystal ordering. This shape-dynamics 
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(a) (b) (c) 



(d) 

Figure 4. Results of growth of an initially homogeneous nematic spherulite in an isotropic matrix in three- 
dimensional simulation where the isosurface corresponds to S = 0.18 and the surface shading corresponds 
to the biaxial order parameter P (minimum/maximum 0/0.07, see colour bar): (a) y/z-plane (b) x/y-plane 
(c) x/z-plane (d) colour bar for (a)-(c) ; phenomenological parameters were n = 0.084 N J^ 2 S , T/vi = 307. 2K, 
a = 1.4 x 10 5 J/m 3 -K, b = 1.8 x 10 7 J/ra 3 , c = 3.6 x 10 6 J/ra 3 , h = 3.0 x lO" 12 J/m, h = 3.1 x lO" 12 J/ra, 
I3 = 0.0 x 10 _ 12 J/ra (pjq) : horizontal length scale is 225nra, and a C++ finite element library, LibMesh 
0.63 was used to develop the fully adaptive parallelised finite-element code. 



regime of growth transitions from the initial nucleus texture to a spherulite with 
shape and texture that minimises the overall free energy, allowing for constant 
growth to proceed. Once this transition is complete, a self-similar growth regime 
is observed where spherulite texture and shape are independent of length scale 
and spherulite dimension scales with time r oc t n . For the liquid crystal 5CB 
(pentyl-cyanobiphenyl), which exhibits an isotropic/nematic transition, the self- 
similar growth regimes for initially homogeneous spherulites were found to range 
from time-scales of approximately and length scales of approximately 1.5/jm 
depending on temperature quench depth (0, 28). 



Past work focused on an initially radially oriented nucleus (1 1% ) , where the initial 
nucleus shape has homogeneous interface conditions and the self-similar growth 
transition involves spherulite core dynamics only. A homogeneously oriented nu- 
cleus has interfacial heterogeneities which result in bulk texture dynamics min- 
imising total free energy, depending on the existence of preferred anchoring at the 
isotropic/smectic-A interface. This results in a prolonged shape-dynamic growth 
regime which has been found to involve defect shedding in nematic cases (li?a ). The 
complex smectic "batonnet" structures composed of curvature defects first observed 



by Friedel and Grandjean (I3il43<5h have been attributed to growth from initially 
homogeneous smectic- A nuclei Past approaches to studying these types of 

growth processes have involved highly simplified shape equations that take into ac- 
count approximations of anisotropic interfacial anchoring energy and bulk energy 
contributions. The simplest approach is that of the Wulf-construction ([H) which 
determines surface shape by minimising the sum of total interfacial energy and an 
ideal undistorted bulk contribution: 



/ -y(r)dA + aV (10) 

J A 



where F is the total free energy of the spherulite, 7 is the interfacial tension (a 
function of position r), and a the free energy density of the spherulite bulk. While 
this approximation is suitable for crystal growth where a homogeneous bulk texture 
is a valid assumption, in the case of liquid crystals the occurrence of bulk elastic 
distortions and complexity of interfacial anchoring energies requires more rigorous 
approaches. However, the Wulf approximation is a convenient starting point for 
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Figure 5. (a) Schematic of the Wulf construction result for the shape of a growing smectic-A spherulite 
with no preferred interfacial anchoring and an initially homogeneous bulk texture (b) log-log plot of the 
spherulite (major axis) radius versus time for the isotropic nematic elastic case (c) power-law fit of the 
spherulite (major axis) radius for the isotropic nematic elastic case. 



analysis of the first simulation case with isotropic nematic elasticity. In the absence 
of smectic-A ordering, the nematic contribution to the interfacial energy, neglecting 
curvature and biaxiality, is (|3d ): 

6V3/i+/ 2 /2 + 3/ 2 (n>fc) 2 /2 

7 486c 5 / 3 1 j 

where k is the unit vector normal to the interface. Equation [TT] shows that in the 
case of nematic elastic isotropy (I2 = Z3 = 0), the nematic contribution to the 
interfacial tension is approximately isotropic and there is no preferred interfacial 
anchoring of the nematic director. Based upon this premise, the Wulf construction 
prediction for this situation should be valid in that interface anchoring conditions 
do not frustrate the bulk. Figure [5fct shows the layer morphology of this scenario 
where a spherulite grows with a homogeneous texture and perfectly circular scale- 
independent shape. 

Figures [Sb-c show simulation results of the long axis and its power law fit versus 
time for the isotropic nematic elastic case. These results indicate that a transition 
from shape-dynamic to self-similar growth occurs on the order of lfis after nucle- 
ation. These results are in agreement with past work on initially textured smectic-A 
spherulite growth (1 12 ) where in both cases the bulk elasticity is dominant over in- 
terfacial energy. As previously mentioned, similar studies of nematic spherulite 
growth show that the transition from shape-dynamic to self-similar regimes occurs 
on a time scale on the order of 80 jis (37) due to the relatively lower magnitude 
of bulk elasticity of the nematic phase compared to interfacial anchoring energy. 
An important conclusion from past work on nematic spherulite growth, that mul- 
tiple shape-dynamic/self-similar regimes occur pre/post-shedding of defects (0), 
which implies that similar phenomena could occur for smectic-A phase-ordering 
at length scales computationally unavailable. In the smectic-A phase, disclination 
shedding events at the isotropic/smectic-A interface in conjunction with bulk tex- 
turing, could result in curvature defects such as focal conic domains observed in 
smectic-A batonnets (L2il). 

Figures EK-c shows the transient morphology of the simulation results for the 
isotropic nematic elastic case which does confirm that there is no substantial pre- 
ferred interfacial anchoring from nematic elasticity (figure EJl): 



fNe = \h(VQ) 2 



(12) 
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(c) 
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time (s) x io" 5 



(d) 



(e) 



(0 
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Figure 6. Simulation results for the isotropic nematic elastic case where the surface corresponds to Re(^) 
(minimum/maximum —0.18/0.18 corresponding to black/white) for simulations times (a) 6.7/is (b) 15.0/is 
(c) 21.0/is (d) surface plot of the nematic elastic contribution (see equation 1121 minimum/maximum 
0.0/2.151 x 10 4 J/m 3 , see colour bar) (e) surface plot of the smectic-A elastic contribution (see equa- 
tion ^3 minimum/maximum 0.0/2.151 x 10 4 J/m 3 ) (f) aspect ratio of y-/x- spherulite axes versus time 
for the isotropic nematic elastic case (g) colour bar for (d)-(e). 



However, the prediction of the Wulf construction (figure^) is not found to be valid 
due to presence texturing in the bulk of the spherulite. An undulation instability 
(fgjl ) is observed along the centreline of the spherulite parallel to the layer normal. 
Figure [6^ shows the smectic-A elastic contributions: 

fAe = \bi |V*| 2 + - A b 2 |V 2 *| 2 - l -eQ : (V*) (V**) (13) 

which indicates that this is a bulk phenomenon due to the gradient of the smectic 
order across the body of the spherulite. As smectic order decreases approaching the 
interface with the isotropic phase, both the free energy penalty for layer dilation 
decreases (due to decreased smectic-A order) and the layer spacing increases (also 
due to decreased smectic-A order, see figure [2]). The relationship was determined 
for a distortion- free smectic-A domain to be 

[2eS A -3&i\~* 

where do is the smectic-A layer spacing, Sa is the nematic scalar order parameter 
in the smectic-A phase. Due to symmetry, the core region of the spherulite has 
maximum smectic order and thus a minimum layer dilation. This overall bulk 
layer dilation results in an undulation instability to dissipate higher energy layer 
expansion via low energy layer curvature §38 ). 

This layer dilation undulation instability, referred to as Helfrich-Huralt undula- 
tions (fgjl . [gj ). has been observed in films of lamellar systems including smectic-A 
and cholesteric liquid crystals under an external layer dilation force §38 ). Figure [3 
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Figure 7. Past simulation results, using the same material parameters in this work (see figure [3j, of 
undulation instabilities in a two-dimensional thin film geometry where the surface corresponds to Re(^). 
The horizontal boundaries correspond to liquid crystal/solid interfaces and the vertical boundaries are 
periodic; the length scale of the domain is approximate 108nra with an initially imposed undeformed 
texture of 25 smectic-A layers (4nm equilibrium layer spacing); adapted from ref. JTJ) 



shows a representative example of past results for the currently simulated model 
system in a two-dimensional thin film geometry with layer dilation imposed by 
perturbing the plate separation distance. As derived using linearised lamellar elas- 
ticity (?), the optimal wavelength for the undulation distortion in the thin film 
geometry (figure [Tj) is approximately (fg3 ): 

Xu = 2n l^-j (15) 

where A is the characteristic length of the lamellar ordering (on the order of the 
layer spacing), and Xu is the characteristic length of the undulation instability, 
and L is the externally imposed length scale. Numerical simulation results shown 
in figure [6k and the ideal layer spacing dependence equation fT4l shows the source of 
the layer dilation driving the instability. Past work studying a quasi-lamellar sys- 
tem, the cholesteric mesophase, has shown that cholesteric pitch gradients (similar 
to the smectic-A layer spacing) are imposed by the presence of the interface with 
the isotropic phase for this system as well (iC). These past results along with the 
current observations imply that a similar undulation instability should be seen in 
growing cholesteric spherulites, although experimental evidence of this phenomena 
has not been observed. This can be explained based upon comparing magnitudes 
of Xu and spherulite radii for both cases, using estimates of their characteristic 
lengths. For typical smectic-A liquid crystals, A is on the order of nm which results 
in Xu < R in the self-similar regime. Thus an undulation instability is predicted 
to be observed in this growth scenario as is confirmed with the present simula- 
tions. For cholesteric liquid crystals, the characteristic length is on the order of fim 
which is also on the order of the maximum spherulite size that has been observed 
experimentally. Thus Xu > L which explains why the growth-induced undulation 
instability observed in the smectic-A system is not observed in the quasi-lamellar 
cholesteric liquid crystal. 

An additional effect of this texture induced by the introduction of smectic-A 
order is that the spherulite shape deviates from the Wulf construction prediction 
of perfectly spherical growth. Instead, the spherulite has a long axis parallel to 
the smectic layers in order to minimise the over spherulite area with high layer 
dilation. The time evolution of the spherulite aspect ratio, shown in figure [6f shows 
a convergence to a value slight below 1 resulting from an equilibrium bulk texture 
without the presence of frustration from interfacial anchoring, but instead from the 
initial homogeneous texture and order parameter gradients induced by the presence 
of the interface. 

For the equal splay-bend nematic case, equation fTTI shows that the nematic con- 
tribution to the interfacial anchoring prefers a planar orientation to the isotropic 
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Figure 8. (a) spherulite shape/morphology predicted by Fournier and Durand's extension to the Wulf 
construction (based on figure 6b from ref. (34)) (b) log-log plot of spherulite (major axis) radius versus 
time for the isotropic nematic elasticity case (c) power-law fit of spherulite (major axis) radius versus for 
the isotropic nematic elasticity case 



interface (average orientational axis parallel to the interface). This scenario is in 
agreement with experimental observations of preferred planar anchoring for the 
isotropic/smectic-A interface (0). In this case the Wulf-construction is not ade- 
quate in that frustration of the interfacial anchoring with the bulk texture will re- 
sult in deviation from homogeneous bulk ordering. A simple extension of Wulf 's ap- 
proach was made by Fournier and Durand which takes into account finite smectic-A 
elasticity \3$ \\ 

F = [ -ydA+ [ f d dV + aV (16) 

J A JV 

As with the Wulf-construction, minimisation of the total free energy of the 
spherulite F predicts its shape, but an additional term is used to describe all 
possible discrepancies with respect to the ideal homogeneous bulk texture: 

f d (r) = f el (r) + SfW(r-A q )} (17) 

where f e \ is the elastic free energy density (curvature/dilation) and the second term 
describe melting of smectic-A order due to defects. Fournier and Durand determine 
a semi-quantitative approximate solution to equations [TBlfTTl which predicts a re- 
laxed configuration as shown in figure This solution determines a spherulite 
configuration where, to a first approximation, dilation, curvature, and interfacial 
anchoring energies are minimised. Non-dilative configurations exist which involve 
complex curvature defects, or focal conies, in order to minimise total energy of 
the growing spherulite (0). Fournier and Durand determine a specific spherulite 
focal conic texture which results in a total free energy substantially lower than that 
predicted by the extended Wulf's method (0). This type of approach is limited in 
that predictions can be made, but transition mechanisms and alternate dilation- 
minimum modes are not energetically identified from a growth process. Thus pure 
geometrical approaches, or those extended to account for dilation/curvature elas- 
ticity are not able to take into account an adequate subset of the physics involved 
in these growth processes. 

Figures 0>c show simulation results of the long axis and power law fit versus 
time for the nematic elastic isotropy case superimposed on those results from the 
nematic elastic isotropy case. The introduction of anisotropic interfacial tension 
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(d) (e) (f) 



(g) 

Figure 9. Simulation results for the equal nematic splay-bend case where the surface corresponds to Re(^) 
(minimum/maximum —0.18/0.18 corresponding to black/white) for simulations times (a) 6.4/is (b) 13.7fis 
(c) 19.6/is. (d) surface plot of the nematic elastic contribution (see equation (|18|) . minimum/maximum 
0.0/2.151 x 10 4 J/ra 3 , see colour bar) (e) surface plot of the smectic-A elastic contribution (see equation 
(|13|) . minimum/maximum 0.0/2.151 X 10 4 J/m 3 ) (f) aspect ratio of y-/x spherulite axes versus time for 
both simulation cases (g) colour bar for (d)-(e). 



due to the equal bend-splay elasticity (figure [6]i): 

ZiVe = ^i(VQ) 2 + ^ 2 (V-Q) 2 (18) 

results in little change in the growth kinetics compared to the previous case. Figures 
EK-c show the dynamic morphology of the spherulite which is in good agreement 
with the shape prediction of Fournier and Durand (figure [Sfci) . Again, the undula- 
tion instability is present in addition to interfacial heterogeneity from the preferred 
planar anchoring. The additional smectic layer curvature resulting from the tran- 
sition from homeotropic to planar anchoring (molecular axis perpendicular to the 
interface) at the spherulite poles also promotes the undulation instability resulting 
in an increased amplitude compared a spherulite of equal vertical radius in the 
previous nematic elastic case. 

Figure shows that the evolution of the aspect ratio differs compared to the 
previous nematic elastic case, where now the aspect ratio decays as the spherulite 
grows from a maximum value following the shape-dynamic growth regime. While 
the spherulite shape morphology and preferred anchoring of the isotropic/smectic- 
A interface are in agreement with past theoretical and experimental observations, 
respectively, the decay of the aspect ratio is not in-line with experimental obser- 
vation of high-aspect ratio batonnets. This decay while in the self-similar growth 
regime can be explained based upon a scaling theory derived from the studying 



shape and director-field transformation of nematic tactoids (l^il Uffl - This approach 
uses a simplified nematic spherulite free energy which takes into account both in- 
terfacial energy coupled with bulk elastic anchoring as a function of the nematic 
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director: 



F = t J (l + u;(q.n) 2 )dA+^K J ((V • nf + [n x (V x n)] 2 ) dV (19) 

where r is interfacial tension, uj is a dimensionless anchoring strength, g is the 
surface normal (a function of position r), n is the nematic director, and K is 
the elastic constant (using an equal splay-bend assumption). Assuming a fixed 
bipolar nematic texture (see figure where contour lines/smectic layers indicate 
the orientation of the vector perpendicular to the nematic director n) and a fixed 
spherulite volume V, a scaling estimate for the aspect ratio of the spherulite obeys 
the following relationship to minimise total free energy: 

^ « ^3/5 r -3/5 y -l/5 (2()) 

r 

where R is the major axis and r is the minor axis. Equation ([201) predicts that the 
spherulite aspect ratio decreases with volume, as is found for the equal splay-bend 
case in figure EJ. This implies that as the spherulite radius grows the aspect ratio 
will converge to 1, assuming no additional shape-dynamic events, which is not in 
agreement with observations of complex batonnet structures in isotropic/smectic- 
A transitions. Thus it is expected that if computational resources allowed access to 
larger spherulite length scales, another shape-dynamics regime would be observed 
similar to defect shedding events seen during the growth of nematic spherulites 



5. Conclusions 

The study of the two-dimensional growth of an initially homogeneous smectic-A 
spherulite in an isotropic matrix was performed via modeling and simulation. A 
mesoscale Landau-de Gennes type model was used that takes into account the cou- 
pling between orientational (nematic) and one-dimensional translational (smectic- 
A) order where layer spacing and the coupling of the nematic director and wave- 
vector are energetically imposed. The effect of two different nematic elasticity con- 
ditions corresponding to isotropic and preferred planar interfacial anchoring were 
performed and compared to past experimental and theoretical studies: 

• A growth-induced dilative layer undulation instability was observed unique 
to initially homogeneous smectic-A spherulite growth. This instability was shown 
to be independent of interfacial anchoring and a result of gradients in the bulk 
order. 

• The use of nematic equal bend-splay elastic conditions was shown to be re- 
quired to adequately model the preferred planar interfacial anchoring experi- 
mentally observed for the isotropic/smectic-A interface. 

• The use of nematic equal bend- splay elastic conditions was shown to result 
in a spherulite shape in agreement with past theoretical predictions for both 
smectic-A spherulites and similarly textured nematic spherulites. 

• An aspect ratio decay of the spherulite was observed in the self-similar growth 
regime, which implies that a second shape-dynamics regime follows at greater 
spherulite radius in order to agree with experimental observations of smectic- 
A batonnets with high aspect ratio. A defect shedding shape-dynamics process 
is proposed, similar to that seen in growing nematic spherulites (0??L 0), to 
transition the bipolar oriented smectic-A spherulite to one in the diverse set of 
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batonnet structures composed of focal conic curvature defects. 

This and past simulation work using the high-order Landau-de Gennes 

type phenomenological model of Mukherjee, Pleiner, and Brand U Q) has shown 
great promise for studying both the isotropic/smectic-A transition and the over- 
all smectic-A mesophase. The key aspects of this model, in addition to its phe- 
nomenological nature, are that energetic couplings are employed between ori- 
entational/translational (nematic/smectic-A) order including: bulk nematic and 
smectic-A order (S/if)), average molecular axis and smectic layer normal (n/a), 
and bulk nematic/smectic-A order and layer spacing (S / if) / \a\) . These couplings 
and the use of full tensorial and complex order parameters capture a sub-set of 
the fundamental physics involved in the isotropic/smectic-A transition unavailable 
using more simplified approaches. 

The main limitation, as with most modeling approaches, lies in the limits imposed 
computationally. Cutting-edge scientific computing approaches such as fully adap- 
tive parallelised finite element libraries (for example, LibMesh (fgfl ) and PETSc- 
FEM (Uj)) maximize the impact of currently available computational resources. 
Through the use of these existing numerical approaches, development of paral- 
lelised post-processing code, and the presently used high-order model, smectic-A 
phenomena at macroscopic lengths scales could be accessed to study diverse multi- 
scale growth phenomena such as the fascinating batonnet structures of Friedel and 



Grandjean iSiASS). 
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